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^ Abstract: In 2009 Hofava proposed a power-counting renornializable quantum gravity 

'— ' theory. Afterwards a term in the action that softly violates the detailed balance condition 

fSJ has been considered with the attempt of obtaining a more realistic theory in its IR- limit. 

^ This term is proportional to uw-^' , where w is a constant parameter and W^' is the spatial 

■^ Ricci scalar. In this paper we derive constraints on this IR-modified Hofava theory using 

J;^ the late-time cosmic accelerating expansion observations. We obtain a lower bound of |ct;| 
that is nontrivial and depends on Avi^, the cosmological constant of the three dimensional 

^) spatial action in the Hofava gravity. We find that to preserve the detailed balance condition, 

O one needs to fine-tune kw such that -2.29 x 10"'' < {c^ ky/) / {HIVLd^) - 2 < 0, where Hq 

• • and r^DE are the Hubble parameter and dark energy density fraction in the present epoch, 

• ^H respectively. On the other hand, if we do not insist on the detailed balance condition, then 

rS the valid region for kw is much relaxed to -0.39 < {c^ kw) / {H19.-q^) - 2 < 0.12. We find 

C^ that although the detailed balance condition cannot be ruled out, it is strongly disfavored. 
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1. Introduction 

In 2009 Hofava proposed a power-counting renormalizable quantum gravity theory, which 
is non-relativistic in the high-energy, or UV, hmit and recovers the Lorentz symmetry in the 
low-energy, or IR, limit [1, 2]. Various aspects of this theory have been widely pursued [3- 
30]. In addition to the theory itself, its various implications in cosmology have also drawn 
wide attention [31-34]. Specifically, there has been a large amount of effort in studying 
the cosmological perturbations [35-43], black hole properties [44-54], gravitational waves 
[55-59], dark energy phenomenology [60-63], observational constraints on the parameters 
of the theory [64, 65], and many others. 

With regard to the background space-time expansion, Hofava gravity in principle mod- 
ifies the GR Friedmann equation with additional terms stemming from its non-conventional 
gravity, thereby contributing to the dark sector. Accordingly, the current and future cos- 
mic observations may provide significant constraints on Hofava gravity, especially when 
connecting it with the cosmic accelerating expansion. 

In this context it was found that the effective Friedmann equation derived from Hofava 
gravity with the detailed balance condition would include a cosmological constant (CC) 
term and a radiation-like a~^ term [31, 32], where a is the scale factor in the Friedmann- 
Robertson- Walker (FRW) metric. This radiation-like term originates from the fourth-order 
spatial derivative terms in the Hofava action. On the other hand, if the detailed balance 
condition is violated, then there would be other dark terms induced. 

Although the gravity action of the Hofava theory with the detailed balance condition 
recovers that of GR in the IR limit, its solution may not be so. For example, in [66] it 
was shown that the black hole solution does not recover the usual AdS-Schwarzschild black 



hole solution in GR. In order to attain a more desirable IR behavior without abandoning 
the simplicity provided by the detailed balance condition in the UV limit, several authors 
introduced a soft violation of the detailed balance condition [1, 67, 68], namely, a term 
proportional to the spatial three-curvature R'^^'. In particular, an IR-modified Hofava the- 
ory that accommodates flat Minkowski vacuum was studied in [68] . The exact solutions of 
spherical symmetry with and without matter were obtained in [69]. 

A question then naturally arises: to what extent can the Hofava gravity violate the 
detailed balance condition? One possible means to address this question would be to 
derive constraints on the IR-modification terms and those obeying the detailed balance 
condition from cosmological observations. For this purpose we consider a cosmological 
model studied in [60] based on the IR-modified Hofava gravity, with the FRW metric 
describing the background space-time and with the energy content that includes radiation 
and dust matter. 

In [60] Park showed that the Friedmann equation of this cosmological model contains 
additional a~^, a~^, and CC terms beyond that in GR. Park identified these terms as the 
effective dark energy (DE) that is responsible for the cosmic acceleration. The observa- 
tions about the expansion history can in principle constrain the behavior of the (effective) 
dark energy and thereby constrain the Hofava gravity. To constrain this Hofava Effective 
Dark Energy {HEDE) by observations, an efficient approach is via a phenomenological 
parametrization of the relevant physical quantities that have been well studied. Once the 
relation between the model parameters and the phenomenological parameters is estab- 
lished, the constraints on the model can be obtained from those on the phenomenological 
parameters that have been derived from observations. Park considered the widely used 
Chevallier-Polarski-Linder (GPL) parameterization of the equation of state of dark energy 
[70, 71], 

WdeIo) =Pde/Pde = tfO + lfa(l -a), (1-1) 

where the constraints on the phenomenological parameters wq and Wa from the updated 
observations have been well studied (see, for example, [72-76]). 

Park [60] explored the feasibility of HEDE by considering three best-fit values of 
{wq^Wo) obtained in [72, 73] where a non-flat universe was considered. It was suggested 
by Park that the existence of some HEDE models that satisfy these three best-fit values 
indicates the validity of HEDE. We note, however, that in principle the dark energy den- 
sity of HEDE is determined once the values of {wq^Wo) are given. As will be shown in the 
present paper, these three best-fit models predict the dark energy densities that are much 
smaller than that required by observations, and have thus already been ruled out. 

In the present paper we pursue a more comprehensive test of the IR-modified Hofava 
theory and its resultant HEDE based on the current observations. We particularly em- 
phasize that a complete test of HEDE based on the cosmic expansion must take into 
consideration not only the evolution of the dark energy density /5de(o)) which involves both 
the present value of Ode and the equation of state Wde(o)5 but also the present value of 
the fractional density ^^ of the spatial curvature. 



Specifically, the present density fraction of dark energy JIde should be around 0.74 
and that of the spatial curvature ^k < 0.01 [75]. In addition, since in the dark sector the 
radiation-like term would be dominant in the early universe, its energy density must be 
smaller than the true radiation energy density $7^, otherwise HEDE would predict a later 
epoch of the matter-radiation equality and that in turn would be in conflict with the cos- 
mic microwave background (CMB) and the big-bang nucleosynthesis (BBN) observational 
results. As will be shown in the present paper, the observational constraint on the dark 
energy equation of state Ti^oE, together with the above three observational requirements, 
very tightly constrains HEDE. These observational constraints on HEDE favor the viola- 
tion of detailed balance of Hofava gravity. We will present a lower bound on the extent of 
the violation. 

This paper is organized as follows. In Sec. 2 we give a brief review of HEDE, a cos- 
mological model based on the IR- modified Hofava theory. In Sec. 3 we discuss the general 
strategy that we utilize for the model test, which involves an approximate relation between 
the model parameters and the phenomenological parameters. In Sec. 4 we investigate the 
observational constraints on HEDE that are presented in the phenomenological parameter 
space. In Sec. 5 we show how the constraints on the phenomenological parameters are tran- 
scribed into that on the model parameters, and investigate its impact on the IR-modified 
Hofava gravity. In Sec. 6 we show the evolution patterns of the effective dark energy in 
the HEDE models that are consistent with observations. We conclude in Sec. 7. 

2. Setup of the Model 

To be self-contained, in this section we give a brief review of the IR-modified Hof ava-gravity 
cosmological model that was investigated by Park [60]. The action of the IR-modified 
Hofava gravity reads 

S, = I dtcCx^jN \^^ (K„K" - \K') - £,C.,C-' + ||£'*4" V,fl»)'j (2.1) 

where the extrinsic curvature 

Kij = ^ i9i, - V.iVj - V,iV,) (2.2) 

(the dot denotes the time derivative), the Cotton tensor 

^ij ^ ^iki^^ fj^{3)j^ _ li?(3)<5A ^ (2.3) 

and K, A, z^, n, A\y,uj are constant parameters. Note that on the right-hand side of Eq.(2.1) 
the last term proportional to uR^"^' induces the soft violation of the detailed balance con- 
dition. For a homogeneous and isotropic universe we consider a FRW metric of the form 



ds^ = -c^dt^ + a\i) 



1 - kr'^/Rl 



(2.4) 



where k = +1, 0, —1 corresponds to a closed, a flat, and an open universe, respectively, and 
Rq is the radius of spatial curvature of the universe in the present epoch. Assuming that 
the matter contribution is in the form of an ideal fluid with energy density p and pressure 
p, Park obtained [45] 



a \ K? 



a) 6(3A-1) 

,.2 



3kV ( -k^ , 2k{Kw-oj) ,2 



a 



a 6(3A-1) 



8(3A-1) \Rla^ Rla^ 

2^^ + '^)^8(3A^V^-M 



(2.5) 
(2.6) 



where the analytic continuation p? — )• —ji'^ for Kw has been employed [45, 66, 77]. The 
upper (lower) sign corresponds to the case where Aw < (Ayi/ > 0). 

Comparing them with the Einstein equations derived from GR with the FRW metric, 

aV S-kG , , c^k ,^^, 

a A.'kG 

- = --^[(/'m + PDE) + 3(pm+PDE)], (2.8) 

one can connect the Hofava parameters k. A, ^u, A^y, a; with the speed of light c, the Newton's 
constant G, and the effective dark energy density poE and pressure Pde, although the 
connection is not unique. Park defined the fundamental constants c and G as 

2_^^V^M ^_ ^V 

8(3A-1)2''" 167r(3A-l)- ^^^ 

To ensure the positivity of the dark energy density as required by observations, we consider 
the case where A-w > and identified the dark energy density and pressure as 



3c I Hq^I^ 2HqUjQ.j^ 2 

^ 4^ 2 2 

3c4 fH^^l 2H^ujni 



p^^ = T,^^^ -T^--^TP' + ^'wl (2-10) 
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where fi/; = — (c^ k) / (RqHq) . The equation of state parameter is then 

_ Pde _ H^^l + 2c^H^ujQka^ - 3c^A^a^ 

^DE — „ —0 7^402 ^.27^2.0 _2 , o„/1a2 „4- l^-J-^j 



Pde ^H^^l - 6c^H^unka^ + Sc^Alra"^ ' 



From Eq. (2.10) we obtain 



= Oia^^ + Osa-^ + Og, (2.14) 

where the critical energy density pc = {3HQC^)/{8irG) and 



H'o^l ^ n c - ^^^ n _ ^'^ 



vy 



"■ = 2A^2»- "^ = -A^- "» = W""' <'•'" 



As shown by the above formulae, this effective dark energy consists of three components 
that are radiation-hke, curvature-like and CC-like, respectively. It involves three model 
parameters: {$7^,0;, Ayi/}. The first parameter stems from the FRW metric ansatz, and the 
last two from the Hofava gravity action. Note that for a flat universe this effective dark 
energy behaves as a CC and accordingly this model is the same as the CC dark energy 
model in GR, i.e. flat ACDM, which is consistent with all the current observational results. 
In the present paper we will consider a nonzero Qk for possible deviation from ACDM. 

3. Model Parameters and Phenomenological Parameters 

To compare the Hofava Effective Dark Energy (HEDE) model with observational results, 
one effective approach is to employ a phenomenological parametrization of the relevant 
quantities, whose observational constraints have been well studied, as a mediator to facil- 
itate the comparison. The observational constraints on the model can be obtained from 
those on the phenomenologically parametrized quantities by invoking an approximate re- 
lation between the model parameters and the phenomenological parameters. 

In HEDE the model parameters include {ilfc, w, Avi/}, which determine the dark energy 
behavior. That is, the dark energy density pj^^{a), as well as w^Eio-), is a function of 
{Qfc, w, A14/}. Phenomenologically, the evolution of dark energy is determined by the present 
density fraction Qbe and its equation of state parameter Woe- Here we invoke the widely 
used CPL parameterization [70, 71] of w^e in Eq. (1.1): Woe = wq + Wa{l — a). In 
addition, since the spatial curvature il.^ is involved in the model parameters, it should also 
be included in the phenomenological parameter space when connecting to the model space. 
In summary, the phenomenological parameters are {ilkj^DEjWQjWa}- Accordingly, we 
have a three-dimensional model parameter space and a four- dimensional phenomenological 
parameter space. A mapping between them is required for constraining the model via the 
observational constraints on the phenomenological parameters. 

The relation between r^DE and {^k, '^j ^w} can be obtained from Eq. (2.13) by setting 
a= 1: 



2c2Avy Aw 2Hq 



nBE{^k,to, Aw) = Tffr^ - ^ + ^Aw. (3.1) 



Following Park [60], we connect the model parameters {cj, Aw} with the phenomenological 
parameters by firstly expanding Woe in (2.12) around a = 1 as 

Woe = WoE\a=l " •»^deL=i (l - a) + ■ ■ ■ , (3.2) 

where the prime denotes the derivative with respect to a, and then identify WoE\a=i and 
~'"^deU=i as wq and Wa in the CPL parameterization. As a result, the approximate relation 
between {w,Ai4/} and {wQyWa} reads 

.,(0 :,^ ^ _ H^^l + 2c'H',u;nk - Sc^A^ 

Sc^ink {H^u;nl - 2^Hl^uA^^ + c^uA^^) 

Itiali'fc, W, Aw ) = 7, , (3.41 

2,{H^,nl-2cmlu:nk + c^Al,f 



or equivalently, 



(1 - 2wo - 3w^ - Wa)H^n 



k 



(1 + 4w;o + Su-g + U'a)c2 ^ ' 

The required mapping between {J7fc,cj, A^y} and {^k^^DE, wo,Wa} is given by Eqs. (3.1), 
(3.3) and (3.4). It maps the 3D model space to a 3D hypersurface in the 4D phenomeno- 
logical parameter space. 

Regarding the observational constraints, as stated in Sec. 1, we assume ^be ~ 0.74, 
r^fc < 0.01 [75], ill < ^r, and the constraint on Woe obtained in [74]. For simplicity, we 
fix JIde = 0.74. This reduces the dimension of the phenomenological parameter space 
from four to three: {QkyWo,Wa}, and that of the model space from three to two, i.e. a 
2-dimensional surface in the 3D space {0^,0;, Av^^}. This 2D surface corresponds to the 
relation fi^. = J7fc(a;, A^y) obtained from Eq. (2.13) with poEia = ^)/pc = 0.74. In this case 
the relations in Eqs. (3.1), (3.3) and (3.4) become Qbe = ^de(w,Avi/), wq = wo{uj,Aw) 
and Wa = Wa{io,Aw). They map the reduced 2D model space {^^(a;, Avi/),a;, Ai^} to a 2D 
surface {Qk{^,Aw),'Wo{uj,Aw)iWa{uj,A\Y)} in the reduced 3D phenomenological param- 
eter space {H.kjWQjWa}. Then, the observational constraint on the model is represented 
by the intersection of this 2D surface and the well-studied, observationally allowed re- 
gions of the 3D phenomenological parameter space. We note that the spatial curvature 
is well constrained. The upper bound of \Qk\ is around 0.01 and could be even smaller 
[75]. Thus, the allowed region in the 3D phenomenological parameter space is very thin 
in the Q^ direction. Accordingly, as a good approximation after imposing the constraint 
rjfc < 0.01, we will simply consider the constraints on the {w(),Wa} plane when taking the 
above-mentioned intersection that presents the valid region of the model. 

4. Constraints on the Phenomenological Parameter Space 

X As commented in [60], for the sake of self-consistency one should require A^ > 0. From 
Eq. (3.6), this in turn requires that 

{wa > -1 - iwo - 3wl Wa>{l- 9u;g)/3} 
or {wa < -1 - iwo - 3wl, Wa<{l- 9w^)/3} . (4.1) 

To simplify the following calculations, we define 

A = l + 4wo + 3w^ + Wa, (4.2) 

B = -1 + 9wl + 3wa. (4.3) 

The self-consistency condition in Eq. (4.1) then reads 

{A>0,B>0} or {A<0,B < 0}. (4.4) 



The valid region on the {wQ,Wa} plane for this self-consistency condition is presented in 
Figure 1 by the shaded area. 
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Figure 1: The valid (shaded) region for the self-consistency condition A^ > 0. 

Rewriting Eqs. (3.5) and (3.6) in terms of A and i?, we have 
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3^ + 5- 
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B fHl 



\n 



k\- 



(4.5) 



Substituting Eq. (4.5) into Eq. (3.1), we obtain 



|Ofc| = sgn(yl)— VAB floE, 



(4.6) 



where sgn{A) denotes the sign of A. For a positive dark energy density, this relation 
requires A>Q, which, as combined with Eq. (4.4), leads to 



{^>0,B>0}, 



(4.7) 



thereby excluding the bottom middle shaded area in Figure 1. 

The current observations suggest f^DE ~ 0.74 and \VLk\ < 0.01 [75]. With Eqs. (4.6) 
and (4.7), these two requirements give a stringent constraint on the parameters A and B 
(i.e. wq and Wa): 

0<AB<(^^] ^9.74x10-^ (4.8) 

VV30.74y 

This constraint largely shrinks the allowed region in the parameter space, which is presented 
by the shaded region in Figure 2. 




Figure 2: The shaded area shows the largely reduced allowed region for the observational 
requirements, Ode ^ 0.74 and \Qk\ < 0.01, that lead to the constraint in Eq. (4.8). 



As indicated in Eq. (2.14), the effective dark energy sector of the (IR-modified) HEDE 
model consists of a radiation- like (J^i), a curvature- like {^2)^ and a constant (r^a) term. 
In the early universe when a is small, the radiation-like term dominates the effective dark 
energy sector. Such a term, if too large, would cause the epoch of the matter-radiation 
equality happened at a time that is later than that suggested by the observational results 
about CMB and BBN. Specifically, this term may contribute to the effective relativistic 
degrees of freedom in the early universe. In the CMB analysis the radiation is usually 
subdivided into two categories: (i) photons and (ii) effective neutrinos (including neutrinos 
and other effective relativistic particles). Accordingly, 






n^ + Q.y = f]^(l + 0.2271Afe//), 



(4.9) 



where 0^ is the present radiation energy density fraction, and N^ff is the number of 
effective neutrino species. The WMAP results suggest that fi^ — 8.47 x 10~^ and N^fj ~ 4 
[75]. Accordingly the effective neutrinos have a comparable contribution to the radiation 
energy density. Regarding the radiation-like term from Hofava gravity as a source of the 
effective neutrinos, we obtain an upper bound of ili given by 0.2271iVejjr2^, which is of 
the same order as Vtj.. This leads us to impose the constraint. 



Q.i<Q.r^ 8.47 X 10^^ 



(4.10) 



Substituting Eqs. (4.5) and (4.6) into the definition of Oi in Eq. (2.15), we obtain 
ill = 3|j4|riDE/8, from which and Eq. (4.7) the above constraint requires 



0< A< 



8* If 



3.05 X 10" 



(4.11) 



This tightly constrained region is presented in Figure 3 by the black area which is so narrow 
that it looks like a black curve. In this figure we also show the constraint given in Eq. (4.8), 



which is presented by the gray area. For wq < —1/3, it largely overlaps with the black 
narrow region. 

In addition, in Figure 3 we show the la (long-dashed contour) and the 2a (dot-dashed 
contour) observational constraints of wq and Wa obtained in [74] from the combined data set 
that includes the SN-Ia data from the Constitution Set, the CMB measurement from the 
five-year WMAP, and the BAO measurement from SDSS and 2dFGRS. The intersection of 
all the above-mentioned allowed regions gives the valid IR-modified Hof ava Effective Dark 
Energy model, which is the black narrow region enclosed by the long-dashed (la) or the 
dot-dashed {2a) contour. 




Wo 

Figure 3: Various constraints on the wo-Wa plane for the IR-modified Hofava Effective 
Dark Energy (HEDE) model. The constraints are given in Eqs. (4.8) and (4.11), and in 
[74]. The long-dashed and the dot-dashed contour respectively present the la and the 2a 
constraint on {wQ,Wa}, which is obtained in [74] from the current SN la, CMB and BAO 
data. The narrow gray region presents the constraint in Eq. (4.8) for the requirement 
[ilfcl < 0.01. The narrow black region that nearly overlaps with the narrow gray region 
for wq < —1/3 presents the constraint in Eq. (4.11) for the requirement Oi < i^r- The 
intersection of these three kinds of allowed regions gives the valid HEDE model. The three 
black dots denote the best-fit models Park considered in [60] , none of which is in the valid 
region. 



5. Constraining the IR-modified Hofava Gravity 



To constrain the IR-niodified Hofava gravity, here we transfer the observational constraints 
on the phenomenological parameters to the HEDE model via the mapping between the 
phenomenological parameter space and the model space. This mapping is given in Eqs. 
(3.5) and (3.6), i.e., in Eq. (4.5), where \^k\ is a function of A and B, as given in Eq. (4.6), 
after we fix J^de = 0.74. With Eqs. (4.6) and (4.7) substituted into Eq. (4.5), this mapping 
can be rewritten as 



CO 



a/3 



24 



IB 
A' 



A 



1 



w 



B, 



(5.1) 
(5.2) 



where the two dimensionless parameters w and Aw are defined as 



C^UJ 



UJ 



Hq^DE ' 



A 



c^A 



w 



w 



HqVLd^ 



(5.3) 



and Sfc denotes the sign of 0^. 

In the valid region, i.e. the black narrow region enclosed by the 2a contour in Figure 3, 
we have the constraints < ^4 < 3.05 x 10~^, as required in Eq. (4.10), and 6.44 < B < 8.48. 
From the above mapping and the constraints on A and B, we obtain 



LJ (e, A 



w 



2 



A 



w 



2 +e 




where 



< e= -^<2.29 X lO""^, 
4 



-0.39< AvF-2 < 0.12. 
For more details, the constraint on {a). Aw} can be read as follows. 



(5.4) 

(5.5) 
(5.6) 



-oo < SkCJ < SkOJ{emax) < as Aw < 2 - emax , 
-OO < SkOJ < SkOj{emax) > aS 2 - Emax < Aw < 2 , 

< SkCJ < \Jtmaxl2 as Aw = 2 , 



Aw {Aw — 2) < SfctD < oo as 2 < Aw < 2 + e„ 

< SkCj{emax) < SkOJ < OO aS Aw > 2 + Cmax , 



(5.7) 



where emax = 2.29 x 10 ^. This constraint is presented in Figure 4, where the dark region 
and the light region correspond to s^ = + and s^ = — , i.e. Q^ > and Q^ < 0, respectively. 



10 



lOOr 



50 



*3 




Figure 4: The valid uj-Ky/ region corresponding to the vaUd region in the wo-Wa space. 
The dark and the hght region correspond to fi/j > and rijt < 0, respectively. 

As shown in Figure 4, Ayy is restricted to values between 1.61 and 2.12, and |lD| has a 
nontrivially lower limit but no upper limit. We note that for a) = the allowed region is 
almost a point, specifically. 



-2.29 X 10"^ < A 



w 



2 <0. 



(5.8) 



That is, in the case where detailed balance is preserved (w = 0) we need to fine-tune 
the value of A for the HEDE model to be consistent with observational results. Thus, 
the cosmological test strongly disfavors, although does not rule, the Hofava action that 
preserves detailed balance. On the contrary, for large \Cj\ the full range of A^i/, (1.61, 2.12), 
is allowed. Moreover, there is no upper limit to |a;|. Accordingly, the observational results 
suggest the breaking of the detailed balance condition. 

Note that the curvature-like effective energy term in Eq. (2.14), as originated from the 
soft violation of detailed balance, remains finite when uj goes to infinity. This can be seen 
in the following. 

rifcw Ode , 



^2 



kw 



2,A-B), 



(5.9) 



where Ode ~ 0.74, < A < 3.05 x 10^'' and 6.44 < B < 8.48 for the valid region. We 
emphasize that even though at the action level the magnitude of the soft violation can be 
arbitrarily large with no upper limit, the corresponding effective energy may still be tightly 
constrained. 
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6. Possible Behavior of the Effective Dark Energy 

Recall Eqs. (2.13) and (2.14), the HEDE consists of the radiation-like {Qia~^), the curvature- 
like (ri2a~^) and the constant-like (f^s) sectors: 

= nia^'^ + n2a^'^ + ris. (6.2) 

Substituting Eqs. (3.5), (3.6) and (4.6) into Eq. (2.15), we obtain 

ni = ^^{l + Awo + 3wl + Wa), (6.3) 

o 

n2 = ^^^{-l + 2wo + 3wl + Wa), (6.4) 

n3 = ^{-l + 9wl + 3wa). (6.5) 

o 

Note that rti + il.2 + ^3 = ^de as required. The requirement Ode = 0.74 and the constraint 
on {wo,Wa} give an allowed region in the {01,02)^3} space: a plane Oi + O2 + O3 = 
0.74 bounded by the box {0 < Sli < 8 x 10-^,-0.05 < O2 < 0.15,0.59 < ^3 < 0.78} 
corresponding to the 2a contour on {wojWa}- 

To show the possible evolution patterns of the HEDE, we consider three sample cases, 
DE.l, DE.2 and DE.3, corresponding to three points in the narrow valid region in Figure 
3. The values of {wq, Wa, ^i, ^2, ^3} in these three cases are as follows. 



(6.6) 



The evolution of the energy density in these three cases is shown in Figure 5. These three 
cases have the same density fraction at present, O^e = 0.74, and share similar evolution 
patterns in the late times up to ln(l + z) ~ 0.25. DE.l resembles ACDM at low redshifts. 
Its energy density remains nearly constant for ln(l + z) < 2.0. This is because Qi and ^2 
are both extremely small in this case. In the cases of DE.2 and DE.3 the energy densities 
increase with z rapidly for ln(l + 2;) > 0.25 due to the Q2 term. Generally speaking, it is the 
difference in ^^2 that makes the evolution patterns distinct from each other in the interval 
0.25 < ln(l + z) < 3.5. In the early universe, the radiation-like term in HEDE would 
dominate. Hence for ln(l + z) > 5 the slopes of ln(p/pc) versus ln(l + z) in different cases 
are roughly the same. The value of Oi determines the value of ln{p/pc) at high redshifts. 
The cases with larger Oi have larger energy densities in the early times. Nevertheless, in 
the early times the energy density of HEDE should be smaller than that of radiation, as 
required in Eq. (4.10). 

In Figure 6 we show the redshift dependence of the ratio Pi/pt, where pi stands for 
the energy density of DE.l, DE.2, DE.3, matter or radiation, and pt for the total energy 
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5.54 X 10-2 0.685 


DE.3 


-0.90 


1.70 X 10-1 


3.18 X 10-6 


1.11 X 10-1 0.629 




Figure 5: The energy density evolution patterns in three sample cases. The solid, the dot- 
dashed and the double-dot-dashed line correspond to DE.l, DE.2 and DE.3, respectively. 



Matter 



Radiation 




Figure 6: The redshift dependence of the contribution of different energy contents to the 
total energy density. The solid, dot-dashed, double-dot-dashed, long-dashed, and short- 
dashed lines correspond to DE.l, DE.2, DE.3, matter and radiation, respectively. 

density. The evolution of the three cases nearly coincide with each other after the time 
when ln(l + z) « 0.5 (i.e. z ~ 0.65), around which the crossing between the HEDE and 
the matter energy density at low redshift happens. This is because after the crossing the 
constant-like (ila) term dominates the dark sector, while these three cases have similar Q^. 
Before this crossing, the three cases behave differently. For ln(l + z) < 4 the contribution of 
the HEDE to the total energy decreases with z, and the decreasing is more rapid in DE.l 
than in DE.2 and DE.3. The main reason is that Q2 in DE.l is negative. In general, 
Pve/ Pt decreases with z more rapidly for smaller 02- For ln(l + z) > 4 (i.e., for z > 50), the 
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contribution of the HEDE to the total energy increases with z. This is a reflection of the 
dominance of the radiation- like (Oi) term over the other two terms in the dark sector. The 
slop of this increment and the contribution of the HEDE at high redshifts are determined 
by the value of Vii. In the case with larger ili, poE increases with z more rapidly and p-o^j pt 
is larger for ln(l + z) > 4. However, the contribution of the REDE would be smaller than 
that of radiation in the early times, as required in Eq. (4.10). 
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Figure 7: The redshift dependence of the contribution of the three dark terms in the 
REDE energy density for the three cases, DE.l, DE.2 and DE.3. The solid, the long- 
dashed, and the short-dashed line correspond to the constant-like, the curvature-like, and 
the radiation-like term, respectively. 

In Figure 7 we show the redshift dependence of the contribution of the three terms, 
{ili(l + z)^, ^2(1 + ^)^, ^3}, in the REDE energy density for the three cases in Eq. (6.6). 
For all the three cases, the constant-like (r^s) term dominates the dark sector in the late 
times, while the radiation- like (r^i) term dominates in the early times. Whether there 
exists a period of the dominance of the curvature-like (^2) term depends on the magnitude 
of ^2- For DE.l, the contribution of the term ^^2(1 + z^'^ is much smaller than the other 
two terms. Hence in Figure 7a the line corresponding to this term cannot be seen. For 
the DE.2 and DE.3 cases, the curvature-like term dominates the dark sector during a 
"middle age" . The precise period of such "middle age" varies from case to case. In general. 
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for larger il.2, the dominance of the curvature-hke term starts earher and ends later, with 
larger contribution from this term to the HEDE. 

7. Summary and Discussion 

In this paper we test the IR-modified Hof ava gravity from the cosmological point of view, in 
particular, from the viewpoint of cosmic expansion. We conclude that the Hoi^ava gravity 
with soft violation of the detailed balance condition is consistent with the current observa- 
tional results on the expansion history of the universe. Specifically, this gravity theory can 
generate the late-time cosmic acceleration with the behavior that is well consistent with 
observations. We note that the Hofava gravity with the detailed balance condition, though 
not ruled out, requires fine-tuning Aw such that —2.29 x 10~^ < {c^ Aw) / (Hq^de) — 2 < 
in order to fit the observational data. This result, together with previous studies [66, 78], 
suggests that the breaking of the detailed balance condition, at least softly, is necessary to 
render Hofava gravity a more realistic IR-limit. 

We obtained the observational constraints on two model parameters, Aw and oj, i.e., 
the cosmological constant of the three dimensional Einstein-Hilbert action and the coeffi- 
cient of the soft violation term. The parameter Aw is well-constrained and it should be 
of the order of the inverse square of the Hubble length, Hq/c^. More precisely, we found 
that Aw, in units of (Hq/c^), is bounded within a small range, (1.6ir2DE5 2.12riDE)5 i-e., 
(1.19, 1.57) after imposing JIde = 0.74. On the other hand, we obtained a lower bound, 
but without an upper bound, to |a;| with regard to the extent of the soft violation of the 
detailed balance condition. The lower bound depends on Aw, and in most cases it is also 
around the order of the inverse square of the Hubble length, Hq/c^. 

With our more comprehensive investigation into the cosmic-expansion test of Hofava 
gravity, we found the Hofava effective dark energy (HEDE) much more restrictive than that 
deduced in [60]. Specifically, the allowed parameter space {wq, Wa} is now much smaller. It 
is a narrow strip beside the parabola l+Awo + SwQ+Wa = around {wo,Wa) = (—1, 0). The 
energy density of the HEDE, with different values of the model parameters, can give rise to 
a cosmological constant as well as non-constant behaviors. For the latter cases with non- 
constant energy density, the main difference between the models therein is the evolution 
of the dark energy density in the "middle age" . This is because during the "middle age" 
the curvature- like {^2) term dominates, while the range of allowed ^^2 is not small, namely 
from —0.04 to -1-1.03. As a result the dark energy behavior in the non-constant cases 
almost coincide with each other at high redshifts. This is because in the early times the 
dark energy is dominated by the radiation-like (ili) term that is highly restricted, namely 
ni < 8.47 X 10"^ 

In our analysis we compare HEDE with observational results by invoking the CPL 
parametrization of the equation of state of dark energy: w^e = WQ+Wail — a), as a mediator. 
We then transfer the constraint on the phenomenological parameters {0,^, ^be, 'WQ,Wa} to 
that on the model parameters {^^^,0;, A^} via an approximate relation between {wQ,Wa} 
and {rifcjCj, Avv'}. Naively, since the dimension of the phenomenological parameter space 
is larger than that of the model space, the 4D parameter space seems flexible enough to 



15 



accommodate the mapping from the 3D model space. Nevertheless, to be prudent, the 
validity of this approximation should be checked. In particular, it is important to verify 
that the two energy densities, pnEDE{z; ^k, ^, ^w) and pcpl{z] ^de, wo,Wa), are consistent 
with each other, where phede is given in Eqs. (2.10) or (2.13) and 

PCPL = poe3"'"('^-i)a"3(i+«'o+«'.). (7.1) 

The approximation is valid when the difference between these two energy densities is sig- 
nificantly smaller than the observational accuracy in the relevant redshift range. This 
consistency check requires further investigation. One way to avoid the possible incompat- 
ibility between the model space and the phenomenological parameter space is to use the 
model to fit data directly, e.g., invoking the x^ fitting to obtain the observational con- 
straints on the model parameters {Clk^-^w^^}- This is under our investigation and will be 
reported in our follow-up paper. 

Hof ava gravity, an interesting alternative gravity theory that breaks the Lorentz sym- 
metry, should ideally be constrained by observations and experiments ranging from micro- 
scopic, solar, astronomical, to cosmological scales. From the cosmological point of view, a 
modified gravity theory changes not only the cosmic expansion history but also the struc- 
ture formation. In the present paper we have shown how the IR-modified Hofava gravity 
can be tightly constrained by the observations about the cosmic expansion. In addition 
to the expansion history, we expect the observations about the cosmic structures, such as 
galaxy surveys and weak lensing observations, would also provide important constraints on 
Hoi'ava gravity. This is worthy of further investigations. 
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